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O ■ Abstract 

Q-i We construct time independent configurations of two gravitating elastic 

^ 1 bodies. These configurations either correspond to the two bodies moving in 

bi). a circular orbit around their center of mass or strictly static configurations. 
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1 Introduction 

In the year 1933 Leon Lichtenstein published a book entitled "Gleichgewichtsfig- 
uren rotierender Fliissigkeiten" ("equilibrium figures of rotating fluids") in which 
he collected his work on this topic, published in Mathematische Zeitschrift in the 
years 1905 — 1928. Lichtenstein showed existence of solutions of the nonlinear 
PDEs describing rotating fluids in various configurations under the influence of 
Newtonian gravity. The simplest configurations are: one axisymmetric rigidly ro- 
tating fluid ball and two self-gravitating fluid balls in circular motion around their 
center of gravity. More complicated configurations treated are toroidal rings or 
even families of nested rings. All these problems have in common that the equa- 
tions become time independent in a corotating system and that they express the 
balance of gravitational, centrifugal and pressure forces. Of particular impor- 
tance is that the total force and torque on the bodies must vanish. Lichtensteins 
method is — in modern language — the implicit function theorem in the con- 
text of bifurcation theory. In one class of problems he begins with an exact non 
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rotating solution and finds a nearby solution with small angular velocity. More 
interesting are the cases with several bodies. Starting from two point particles on 
a circular orbit around their center of gravity, one "puts two small homogeneous 
fluid bodies on these orbits": this is no solution to the equation, however for a 
small body the error is small enough to begin an iteration which converges to an 
exact solution. The purpose of this paper is to show that basically all configura- 
tions Lichtenstein treated can also be successfully studied if the bodies are made 
of elastic material. The case of one slowly rotating elastic body was studied in 
[1]. In this paper we treat: Firstly, one elastic body in circular motion in the 
gravitational field of some central mass which is kept fixed. Secondly, two bodies 
in circular motion around their common center of gravity under their mutual 
Newtonian gravitational field. 

We also deal with a static problem, where one small body is near a non-degenerate 
equilibrium point of a given gravitational field - such as that in the center of a 
hollow triaxial cylinder, (iv) like (iii) - with the mutual interaction and self- 
interaction of both bodies taken into account. 

The paper is organized as follows. Section 2 reviews elastostatics; quotes prop- 
erties of the elasticity operator and formulates the equations for general "force 
functionals" . In section 3 we discuss various types of forces, which in particular 
include centifugal and gravitational forces. Section 4 contains the main applica- 
tion, the motion of a body in an external gravitational field and the motion of two 
bodies under their gravitational attraction. Section 5 is in a way separate from 
the main theme of rotating bodies. It shows, that the techniques developed can 
be used to demonstrate the existence of Newtonian static two-body solutions for 
particular geometries of the bodies. This possibility had always been conjectured, 
but the existence of a full PDE solutions had not been demonstrated. 

2 Elastostatics 

In this section we recall the setup of nonrelativistic elastostatics. We are using 
the material representation where elastic states are decribed by deformations, 
i.e. maps from material space, the "body" , to physical space. A fuller and more 
general description can be found in [1]. 

Let Q be a bounded open connected subset in IR 3 with smooth boundary dQ. 
The basic objects are maps $ : Q — > R 3 . We will use X A , resp. x % for Cartesian 
coordinates in Q resp. R 3 and Sab, resp. Sij for the metrics on the domain, resp. 
target space. There is a trivial or reference state, denoted by which we will 
always take to be the identity. We assume the $'s are sufficiently close to $ (see 
below), so that they have a regular inverse / : $ _1 (fi) — > f2. We define by 
the equation 

* A i (X)(d B &)(X) = 6 A B (1) 
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and the tensor H AB by 

H AB {X) = ^ A i (X)^ B J (X)5 ij . (2) 

In other words, H AB is the inverse of the pulled-back-under-$ Euclidean metric 
on R 3 . 

The material is described by the stored-energy function w(X, d$) = w(X, H BC ), 
smooth in its arguments and subject to the conditions 



/ dw \ 

[dHABj 



= (3) 

H CD = S CD 



and that Labcd, defined by 

(4) 



/ d 2 w \ 
Labcd - y dH AB dH CD J 



H EF = S EF 



should satisfy 

Labcd M ab M cd > C\M\ 2 (5) 

with C a positive constant and M AB = M^ AB \ Condition Eq.([3]) will mean 
that the reference state $ has vanishing stress. Condition Eq.fjSJ), called uniform 
pointwise stability, is an ellipticity condition satisfied by standard materials. We 
will sometimes need the condition that the material be homogenous and isotropic. 
Homogeneity means that w is independent of X. Isotropy of the material means 
that w depends on H AB only through the principal invariants of H A b = H ac 5bc- 
In the latter case Labcd is of the form 

Labcd = A Sab^cd + 2/i 5 C (a$b)d, (6) 

and uniform pointwise stability is then equivalent to \i > 0, 3A + 2/i > 0. 
Remark: Perfect fluids can be described by assuming the stored-energy to be a 
function solely of J, where J is the determinant of H A b- However in this case 
condition (jSJ) is never satisfied. 
The "first Piola stress tensor" is defined by 

Let now C be a neighbourhood of $ = 5> in W 2,P (Q, 1R 3 ), p > 3, small enough so 
that each $ e G has a C 1 - inverse (see p. 224 of [4]). 

Next let £ be the quasilinear second-order operator sending each map $ G C to 

£,[*](*) = (W)P0 (8) 
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where ai A {X) = crf(X, d$(X)). Denote <Tj by {(Ji A nA)\aci, where ha is the out- 
ward normal to dVt. Let the "load space" L be defined by L = W°' p (Vl,M 3 ) x 
W 1 ~ 1 ' p ' p (dQ,M 3 ). It is then well known (see [9]), that the operator E sending 
$ G C to the pair {b = Zi^r = at) <E L is well defined and C l . Here are some 
properties of the map E. 

Property 1: E(&) = 0. This is obvious, since <7^(A", <9$) = by Eqs.flSEI)- 
Property 2: If (6, r) G E(Q) then there holds 

/ (f o$)(X) dV(X) + / (f o$)(X) Ti(X) = 0, (9) 

where £*(#) is any Killing vector of Euclidean M 3 . Note that (Q, when required 
only for translations, is independent of $. The conditions ([9]) - often called 
"equilibration conditions" - have the meaning of vanishing total and force and 
total torque. If these conditions were violated, the body will be set into motion 
and no time independent solution exists. 

Finally, the linear map 5E : W 2 ' P (Q,M. 3 ) — > £>, given by the linearization of E at 
$ = $, has 

Property 3: SE is Fredholm with (6-dimensional) kernel given by the Euclidean 
Killing vectors at the reference state, i.e. <5$* of the form 5^ l (X) = (^ J o$)(X). 
Furthermore SE has a range (of codimension 6) given by the subspace of all 
(6, r) satisfying the equilibration conditions at $ = 2>. 
A proof of Property 3 can be found in [9] . 

We now come to the concept of force density in the present context. Let F be a 
"load map", i.e. a C 1 - map from C to L. We are interested in solving equations 
of the form 

E[$] + AF[$] = (10) 
for small A. Our load maps will be of the form 

F:$ee^(^ = ? j [$] )Ti = 0). (11) 

We are thus seeking solutions to the elastic equations with <jj, the normal stress 
at the boundary, vanishing, i.e. "freely floating bodies" . 

We will call a force field 3 equilibrated at $, if satisfies Eq.flH]) with Tj = 
i.e. if 

/ (Co<t>)(X)^[<t>](X)dV(X) = (12) 
Jn 

for all Euclidean Killing vectors £. Given a functional 9^, we will try to solve 
(flQj) by the implicit function theorem. Suppose $a is a 1-parameter family of 
solutions of (flUj) with $ — ^- Then (1121) must hold for any A, whence for A = 0. 
Therefore the reference configuration must be equilibrated for the force functional 
in order for being the starting solution for an application of the implicit function 
theorem. Suppose we have such a $, then the linearisation of the elasticity 
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operator at $ still has a non trivial finite dimensional kernel and range. In 
the spirit of bifurcation theory one could try to project the equation onto some 
complement of the range. Sometimes this approach works, provided one takes 
"sufficiently small" but finite bodies of any shape. We used this approach in [2]. 
In the present section we will just consider spherical bodies. Then one can use 
some discrete symmetries to make the space of configurations and loads smaller 
so that the implicit function theorem can be applied. 

3 Forces 

In our applications we will deal with the centrifugal force and the gravitational 
force under various circumstances. 

3.1 Spatial force fields 

Both the centrifugal and an external gravitational force give rise to maps £F of 
the form 

? l [$](X) = (K t o$)(X) 1 (13) 
where Ki(x) is a smooth vector field on IR 3 with a potential 

Ki = diK (14) 

satisfying 

AK = 5 ij didj K = const . (15) 

We call such fields harmonic. It follows that where is any Euclidean 
Killing vector, is a harmonic function. Let us define 

mAZ; *)= f i?°*){x) (KiO$){x) dv(x), (16) 

JB R (x) 

where Br(x) is the open ball of radius R centered at the point x. Suppose x is 
a critical point of K. Then, by the (solid form of the) mean-value theorem for 
harmonic functions (see [5]), there holds 

[6 = / ?Ki dV(x) = \B R \ {?Ki)(x ) = (17) 

JB R (x ) 

This means that 3^ = Ki o $ is equilibrated at $ provided that Q = Br(xq). 
We will also use the variational formula 

50/W£;$W$= I iCK\,6v dv {.,■). (is) 

JB r (x ) 

Taking in (118jl for = r/, with rj another Euclidean Killing vector, we find the 
integrand in (TIBl) is again harmonic. Consequently 

<WW£; *]-ri= \Br\ (CrfdidjK^Xo) (19) 
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3.2 Force fields equilibrated at all deformations 

It seems remarkable at first sight, that there should exist such force fields. Sup- 
pose there is a symmetric tensor field 0jj[$](x) on R 3 , having the property that 

djGj = ViGl 3 \ (20) 

and that 0y is of 0(|x| -4 ) at infinity. Then 9^ = J(djQj) o $ is equilibrated 
for all $. (Recall that J denotes the Jacobian determinant of $.) To see this, 
compute 

/ (e°®)(X)J(X)((d J <dS)o$)(X)dV(X)= [ C(x)(d j a i j )(x)dV(x). (21) 
Jn J<s>(n) 

Now the integral on the right-hand side of (I2TT) can be extended to all of R 3 . 
After partial integration, using that £ diverges at most linearly at infinity, there 
follows 

/ Vi[$](X)dV(X) = - [ (d j £ i )(x)G ij (x)dV(x) = (22) 
Jn Jm 3 

An example is afforded by the following situation: Let p be the function on R 3 
given by 

pm(x) = [(Jo^)(x)}- 1 p oXH n ) (23) 

with x be the characteristic function and po a constant. This is the spatial density 
corresponding to the (constant) reference density po- Consider, then, the function 
U on R 3 given by 

AU = —AnGp (24) 

with U(x) = 0(|x| _1 ) at infinity, i.e. the gravitational potential generated by p. 
Now take 

©« = -^[(ftE7)(W - \&<i ^(dkUMU)]- (25) 

Then 

fyQJ = -pdiU (26) 

so that fl20l is satisfied. In this case, equlibration can also be checked by first 
computing 

j mX ) = Gpo jf ^fgrff, *V(n •. = ««•' (27) 

and using the explicit form of the Killing vectors. We have of course rediscovered 
the well-known fact that the total force and the total torque on a body due to 
its own gravitational field is zero. 
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4 Circular orbits 



4.1 Test body in a central gravitational field 



Here we have consider the gravitational field of a point source of mass M, in a 
reference frame rigidly rotating with angular frequence u>. This leads to a spatial, 
harmonic force field given by 



K t = po di 



(x l ) 2 + (x 2 ^ 2 



+ 



GMi 



(21 



where r = a/ (x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 . The field K{ vanishes on the circle given by 
r = L, x 3 = where 

n m 

(29) 



u 2 L 



L 2 



which of course corresponds to circular orbits of point particles. Thus we take 
\3i = XKi o $ with 

-(x l ) 2 + (x 2 ) 2 



Ki = di 



L 3 

r 



(30) 



on R 3 \{0} and solve 



E + XF = 0, 



where E and F are respectively defined after Eq. 
near the reference state defined by 

$ = id|n, 



(31) 

and in Eq. ffTT]) . For small A 

(32) 



where Q = Br(1). Here R < L and I is for concreteness taken to be I = (L, 0, 0). 
(The body is here and henceforth assumed to have spherical shape in its relaxed 
configuration.) The interpretation of a solution $^ fc> r small A will be that of an 
elastic body moving with constant angular frequency u along a circle (of radius 
~ $ 1 (L,0,0)) in the central gravitational field of a mass M. Furthermore there 
holds 



A 



p o L0 



GM 
~1? 



(33) 



Since the function in square brackets in Eq. (l30|) has Laplacian equal to a 2, it 
follows from our previous discussion that 3 is equilibrated at $ = 
Our application of the implicit function theorem will be much simplified if we a 
priori impose the condition that the configurations be symmetric with respect to 
reflection at the planes X 2 = and X 3 = 0. Thus we now restrict W 2 *(Vl;. 
to the Banach space W 2 * (fi,M 3 ) of all maps $ G H/ 2 ' p (fi,M 3 ), so that 



$ 2 (X\X 2 ,X 3 ) 
$ 3 (X\X 2 ,X 3 ) 



^(X\-X 2 ,X 3 ) 
-<f> 2 (X\-X 2 ,X 3 ) 
^ 3 (X\-X 2 ,X 3 ) 



$\X 1 ,X 2 ,-X 3 ) 
$ 2 (X 1 ,X 2 ,-X 3 ) 
-<f> 3 (X 1 iX 2 ,-X 3 ) 



(34) 
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These relations are the same as those satisfied by a vector field, which has mirror 
symmetry w.r. to both the X 2 = and the X 3 = plane. Since the set 
(|34l) contains the identity, there will again be an open neighbourhood of the 
identity, which we call C sym , so that each element of C sym has a C l - inverse. We 
will define as restricted load space the £ sym given by all (f>, r) G W /0,p (fi,lR 3 ) x 
W 1 ~ 1 / p ' p (dfl,, R 3 ), with components having the same symmetries fl34|) . Supposing 
that the stored energy function is homogenous and isotropic, it follows that dA&i A 
also satifies (134"|) . Similarly, using that Q - whence dfl - is invariant under the 
mirror symmetry, it follows that <Ji A riA is also invariant. Thus the elasticity 
operator E maps C sym into £ sym . Next observe that the equilibration conditions 
(Q in C sym are by symmetry all identically satisfied except for that where £ = d±. 
Thus £g ym , the set of equilibrated loads in £ sym , is simply given by all {b, r) G 
£ sym , for which 

f h(X) dV{X) + f n{X) dS{X) = 0. (35) 

JB R (l) JdB R (l) 

In particular the set £>t ym does not depend on $. It is clearly a Banach subspace 
of £ sym of codimension 1. Note next that the only Killing vectors in W^(Q, R 3 ) 
are translations in the X 1 - direction. Thus the linearized operator 8E at $ has 
kernel exactly given by the linear span of 5$ = d\. 

Now turn to operator F given by $ G C sym i— > (6j = o $, n = 0) G £. By well 
known results on "Nemitskii operators" (see [9]), this map is C l . Furthermore it 
restricts to a map from C sym to £ sym . Now consider the subset C C sym subject 
to 

Hm\di;$\ = [ K 1 o<S>d 3 X = 0. (36) 

JB R (l) 

As shown in Sect. (3.1), the identity map $ lies in Cj-. Using ( 1181) and Eq.fl30l). 
an easy computation shows that 

^,/[<9i;<5] -ai = |^l "3 (37) 

Thus Cg- is a C 1 - submanifold C C sy m near $ of codimension 1 (see e.g. [6]). 
Furthermore, by construction, E + XF maps into £sym an( i ^s linearization 
at A = 0, $ = $ is an isomorphism. The implicit function theorem immediately 
gives the following 

Theorem 4.1: Let the stored energy function entering the elasticity map E be 
homogenous and isotropic. Fix a sphere of radius R and a number L > R. Let 
u, A be defined by (|33|) . Then there exists a positive number e - which depends 
on L, R and the stored energy function w - so that a solution $a of Eq. (!3~T|) exists 
for A G [0,e). 



S 



X s 



-x l 



Figure 1: Two spheres in steady rotation 



4.2 Test body with self-gravity 

Now the load jFj is given by 

Po L 3 f ^(X)-^(Y) 



X3; = X 



Ki o $ 



m y n |$(x)-$(y)|3 



dV(Y) 



(31 



with Ki given by Eq. (l3"01) and A given by (1551) . Using Sec(3.2) and the fact that 
(see [2]) Eq. fT2~Tl) defines a C l - map from C sym to <Cg ym C £ sy m, everything literally 
goes through as before and we have the 

Theorem 4.2: Under the same assumptions as in Theorem (4.1) the equation 
E + XF = has a unique solution $a with $o — for sufficiently small A. 



4.3 Two identical bodies with gravitational interaction 
and self-gravity 

Here we take for the relaxed state two spheres of radius R at distance L along 
the X 1 axis. Thus, for the reference configuration we take balls Q = Br(^I) and 
Q' = Br(-^I), where the vector I is again given by I = (L, 0,0) with 2R < L, 
and deformations $ : £1 —> M 3 , $' : tt' —> M 3 with 

^(X 1 ,* 2 ,* 3 ) = -$\-X\X 2 ,X 3 ) 

$ ,2 (X\X 2 ,X 3 ) = <$> 2 (-X\X 2 ,X 3 ) (39) 

$ ,3 (X\X 2 ,X 3 ) = $ 3 (-X\X 2 ,X 3 ) 

Thus the relaxed bodies and their deformations have complete mirror symmetry 
with respect to the plane X 1 = (see Fig.l). 
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The force field is given by 



2|BrI7b s (1„ |*(x)-*(y)|s 

where a = 1,2. The interpretation of a solution $^ is that of a pair of elastic 
bodies of mass m = \Br\ p Q (reduced mass M r ), rotating about their centre of 
mass with angular frequency u where 

o 9 T GM r m 
X = Po lo\ u 2 L=^f, M r = —. (41) 

The self-interaction term in (TJO]) is again equilibrated for all <3>'s. The first two 
terms in (j4T)l) are again of the form Ki o $ with Ki a harmonic force field, where 
however, for the second term, the potential for this force field is in turn a func- 
tional K of namely there holds 

AK[&](x) = -4vrp, (42) 

where (see (123]) ) 

p[$'](x) = [(J'o$'- 1 )(x)]-VoX^ (f , ) (43) 

with J' the Jacobian of the map We again find that 5T is equilibrated at the 
identity. Furthermore 

6* f i Rt i l [d l ;$].d l =\B R \.3 (44) 

for the sum of forces in (1401) . of course the last term does not contribute to /i. 
Thus exactly as in the previous section is a C 1 - submanifold C C sym near 
$ of codimension 1 and E + XF maps into -Csym with its linearization at 
A = 0, $ = $ an isomorphism. We have obtained the 

Theorem 4.3: Consider Eq.( fT0|) with homogenous and isotropic stored energy 
function and with F of the form F : $ 6 C i-> (fej = 9 r i[$],Ti = 0) where 9^ is 
given by Eq.(140l). For sufficiently small A there exists a solution $a with $ = 



5 Static problems 

In Newtonian theory it seems obvious that there exist no static solutions with 
two gravitating bodies: the bodies should fall onto each other. This however is 
in general only true if the geometry is such that the bodies are separated by a 
plane. In this section we show that for more involved geometries there do exist 
static 2-body solutions. 

First we try to put a small body near a critical point of the field of some external 
body. The gravitational force vanishes at critical points of the gravitational 
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(a) (b) 
Figure 2: static configurations 

potential due to this body. In section (5.1) we will simply assume that a critical 
point, subject to some further restrictions, be given, and then show that there 
exists an equilibrium configuration of a small elastic body near that critical point. 
Such critical points are however absent in those parts of the vacuum region which 
are separated by a plane from the support of the large body. One thus has to 
put the small body in the "hollow space" exterior to a sufficiently non-convex 
body (Fig. 2a). We give a particular example for such a body in section (5.2) 
corresponding to Fig. 2b, where we also consider the full 2 body problem. 

5.1 Test body 

We first treat the case of a body in a given potential field Ki = diK. We make 
the following assumptions: 

1. K has a nondegenerate critical point at the origin. 

2. The domain Vt has the origin as centre of mass. 

3. The Hessian of K and the tensor of inertia of Q are simultaneously diago- 
nalizable (i.e. commute). 

4. The eigenvalues of the inertia tensor are all different, and the eigenvalues 
of the Hessian are all different. 

We now show that, by simultaneously translating, rotating and scaling Q, one can 
satisfy the equilibration conditions (at the identity). Let = fuj\ be the Hessian 
of K at the origin, We first observe that the force Ki = fijX^ is equilibrated w.r. 
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to n. Namely 

lin[d l du $] = / d%-xW(x) = 0, (45) 
Jn 

by Ass. 2 and 

fa[u i j xfa i ;$]=u ij f ik e j k , (46) 

where 

Qi k = j x j x k dV{x). (47) 
Jn 

Expression (!4"6l) is zero on account of o;^- = cj^j together with Ass. 3. (The 
quantity O u differs from the standard definition of the inertia tensor by a multiple 
of 5 l i , so properties 3 and 4 carry over to O u .) The full force Ki can be written 

as 

Ki = furf + f ijk (x)x j x k . (48) 

Denote by &s,d the map x l i— > = S^a^ + Z) 1 with 5* a rotation matrix and 
consider the expression 

k{d,e;S,D) = ^vnAd i di;$ S)tD }, (49) 

where f2 e = {x G R 3 |^x G Q}. After the obvious change of integration variable it 
follows that 

k{d,e;S,D) = [ d' f ;j (S{ r! / ■ IP) ■ 
Jn 

+ ef ljk (eSy + eD)(S j f y j ' + D i ){S k k >y k> + D k )}dV(y). (50) 
Similarly, the quantity 

m(u, e; S,D) = /^JuA^'d;; (51) 

is equal to 

m{u,e;S,D) = [ iv i j {S j j ,yi'+Di)[f ik (S k k ,y k ' + D k ) + 
Jn 

+ ef ikl (eSy + eD)(S k k ,y k ' + D k )(S l v y 1 ' + D l )]dV{y). (52) 

Clearly we have that k(d, 0; id, 0) = 0, m(uj, 0; id, 0) = 0. We can view the pair 
(k( ■ , e, S, D), m( ■ , e, S, D)) as defining an e - dependent map from the Euclidean 
group E(3) into the dual of the Lie algebra of E(3), and the derivative of this 



map as a bilinear form V on the Lie algebra of E(3). We find that 

5 D k(d, ; id, 0) • e = \Q\ (Pf^ (53) 

5 s k(d, ; id, 0) • n = (54) 

6 D m{u, ; id, 0) ■ e = (55) 

<y s m(o;,0;id,0)./i = (faG** 1 + f mi <d> ml 5 k j jwV'h (56) 



12 



where /ijj = yU[y] and we have used Ass. 2 in (1541) and (|55|) . The expression in 
brackets on the r.h. side of (|56|) . viewed as a bilinear form on antisymmetric 
tensors, is symmetric by Ass. 3. Writing co l j = e l jkp k , = e l jkq h , this form is 
equivalent to | O | r^p 1 ^ where r l j, by the standard e - identities and using Ass. 3, 
is in matrix notation given by 

T = 3 9 / - (tr /)9 - (tr 6)/ + [(tr 0)(tr /) - 2 tr (6 /)] id (57) 

Going to a frame where both and / are diagonal, we see that r is also diag- 
onal, with e.g. Tn = (/33 — /22)(@22 ~ ©33)- Thus r is non-degenerate, using 
Ass. 4. Using Ass.l, the full bilinear form V is also non-degenerate. There- 
fore by the implicit function theorem we conclude that there exists, for e suffi- 
ciently close to 0, a unique family (S(e), D(e)) with (S'(O), D(0)) = (id, 0) so that 
Hn c [d % di; ®s(e),eD(e)} = for all d % and nn e [u} l jX 3 di] ®s(e),eD(e)} = for all u l j with 

Uij = U)[ij]. 

We next define fl = Q e for fixed small positive e and, by slight abuse of notation, 
5> = idQ. Again our task is to solve (fit)]) , namely 

+ AF[$] = , (58) 

where the load map F is of the form 

F:^ee^(b i = ^],T i = 0) (59) 

and 3 r i[$] = Ki o $. Both E and F are viewed as maps from C to the load space 
C (see Sect. 2). We have the following 

Theorem 5.1: For sufficiently small A equation (1581) has a unique solution $^ £ C 
near $. 

Proof: In a similar way as in Sect. 4.1 we define Qj- as a neighborhood of $ in C 
satisfying the 6 conditions 

ma[£;S] = Af °*)(*) WW ^(*) = o (60) 

for all Euclidean Killing vectors £. Consider the expression 

; = / (£ e F)< 5$ 4 dV(X) (61) 

If (5>$* is any non-zero Killing vector rj and e entering the definition of Q is suf- 
ficiently small, it follows from the above discussion and continuity that the ex- 
pression (IBTp can not be zero for all £. Thus is a C 1 - submanifold of C of 
codimension 6. 

There is an operator P : £ — > projecting onto loads equilibrated at the iden- 
tity. Such an operator is constructed by choosing a (6 - dimensional) complement 
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of £ $ in £ and requiring this complement to be the kernel of P. We now consider 
the modified equations 

Po(E + XF) = 0. (62) 

for $ G e<r. The linearization 5E : W 2 > p (Cl) -> £$ at A = of the map in Eq.(j52j) 
is clearly an isomorphism, so there exists, for small A, a solution $ of Eq. (|62j) 
near $ G C^. By property 2 of the elasticity operator in Sect. 2 we have that, for 
all $ G C, E[<&] G the set of all loads satisfying the equilibration conditions 
at $ = $, so this holds trivially for $. We also have that F[$] G £|,, since 
<£ G ey. Thus (E + AF)[$] G £|>. By Eq.(l62l the load (E + AF)[$] lies in some 
complement of This complement, by continuity, is also a complement of 
Thus (-E 1 + AF)[$] = 0, and the proof is complete. 

5.2 A static 2-body problem 

We first need a domain whose Newtonian potential (assuming constant density) 
has a critical point in the vacuum region with all three eigenvalues of the Hessian 
non-zero and different. 

Consider first the Newtonian potential U of a solid ellipsoid centered at the origin 
with axes a± > a 2 > a 3 and constant density p'. From chapter 3 of [3] (Eq.'s (40) 
and (18)) we read off that that 

(d^UoM = -27rGp , diag(A 1 ,A 2 ,A 3 ), (63) 

where < A\ < A 2 < A3. In the limit when 03 approaches a 2 we find (using either 
Eq.(18) of [3] or the explicit expressions (38,39) of [3]) that A x < A < A 2 = A 3 , 
where A is the value of A4 when a\ = a 2 = 03 (which is equal to |). By continuity 
we have A\ < A < A 2 < A 3 for 03 less than but close to a 2 . We now subtract 
from U the potential of a solid sphere of radius a < a 3 , centered at the origin, 
thus obtaining a new potential U corresponding to a hollow triaxial ellipsoid Q'. 
The potential U has a critical point at the origin with Hessian having one positive 
and two negative eigenvalues, all different. 

We now choose a domain Q with mass center at this critical point and such 
that the inertia tensor satisfier Ass.'s 3 and 4 in Sect. 5.1: as explained in Sect. 
5.1 we can, by translating, rotating and scaling Q, find a domain Q so that 
£*n[£;idn] = for all Euclidean Killing vectors £. By slight abuse of notation 
we henceforth call fl again Q. Using an argument identical to that in section 
(3.2), one sees that when pn[£] id^] = there also follows that pn'[£] idc] = for 
all Euclidean Killing vectors ("actio = reactio") - and this can also be checked 
explicitly. We now assign to the two domains Q' and Q the constant densities 
p' and p and elasticity operators E' mapping $' G C into the load space £/ and 
E, mapping $ G C into the load space £, respectively, with both stored energy 
functions satisfying Eq.'s (131151) . Neither isotropy nor homogeneity of the stored 
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energy functions is needed in this section. Associated with Q' there is the load 
map F' with 3"' given by 

Xm *<](X<) = -G[ P I jgg^fgji dV(Y) + 

+p I i^o-4-)P dV(r) ] (64) 

and with f2 there is associated J given by 

We now define DjCCxC'c W 2 ' P (Q,R 3 ) x VF 2 ' p (f2', R 3 ), p > 3 as the set of 
(<£>, $') G 6 x C for which 

L^ i'w-S ^^' (66) 

and $' satifies 

= 0, d' [l $> ] (X ) = 0, (67) 

where X is some given point in Q'. Again, by section (3.2) or by calculation 
using the explicit form of C(X), one shows that Eq. (1661) is equivalent to 

Combining the fact that there are no non-zero Killing vectors vanishing at a point 
together with their curl and the discussion of the previous subsection - making 
e entering the definition of Q smaller if necessary, it follows that is a C 1 - 
submanifold of C x 6' near ($, $') of codimension 12. We now set A = G and 
solve, using the implicit function theorem, the equations 

F(F[$] + AF[$,$']) = 0, P'(E'[<S>'\ + AF' [$,$']) = (69) 

on Dy, for small G and ($, $') near ($, $'). Here F, F' are projection operators 
on fl resp. £/, defined analogously as in Sect. (5.1). The argument showing that 
we have in fact also solved the "unprojected" equations proceeds as in Sect. (5.1), 
if it is recalled that the self force terms in (I65II64P are identically equilibrated. We 
thus have the 

Theorem 5.2: Consider the coupled set of equations 

F[$] + AF [$,$'] = 0, F'[$'] + AF [$,$']= 0, (70) 

where ($, $') 6 6x6'. Then, for A sufficiently small, there is a unique solution 
($ A ,$' A ) with ($ ,$o) = (idnjidnO- 
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